Theory of parametric amplification in superlattices 
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We consider a high-frequency response of electrons in a single miniband of superlattice subject to 
dc and ac electric fields. We show that Bragg reflections in miniband result in a parametric reso- 
nance which is detectable using ac probe field. We establish theoretical feasibility of phase-sensitive 
THz amplification at the resonance. The parametric amplification does not require operation in 
conditions of negative differential conductance. This prevents a formation of destructive domains 
of high electric field inside the superlattice. 

PACS numbers: 03.65.Sq, 73.50.Mx, 73.21. Cd, 84.30.Le 



Motion of an electronic wavepacket in a periodic lattice 
potential with a period a subject to a constant electric 
field Edc is characterized by oscillations of its velocity 
with the Bloch frequency lob = eaEdc/h [!(■ Bloch oscil- 
lations originate in Bragg reflections of the particle from 
the Brillouin zone boundary. Among solid state struc- 
tures, artificial semiconductor superlattices (SLs) with a 
relatively large period and narrow bands are most suit- 
able for manifestation of Bloch oscillations effects [|J • In 
the stationary transport regime, Bloch oscillations causes 
static negative differential conductivity (NDC) of SL if 
lobt > 1 (t ~ fOO fs is a characteristic scattering time) 
2]. For lubt > 1 and homogeneous distribution of elec- 
tric field inside SL, it canpotentially provide a strong 
gain for THz frequencies [3]. However, in conditions of 
static NDC the same Bragg reflections, which give rise to 
Bloch oscillations, do excite a soft dielectric relaxation 
mode resulting in a formation of domains of high field 
inside SL jij. The electric domains destroy Bloch gain 
in a long SL. Therefore an utilization of Bloch gain is a 
difficult problem 

Using simple semiclassical approach, let us consider 
now an influence of Bragg reflections on dynamics of an 
electron subject to a strong ac (pump) field E p (t) = 
Eq cos tot. Combining the acceleration theorem for the 
electron momentum along the SL axis, p — eE p (t), and 
the tight-binding energy-momentum dispersion for a sin- 
gle miniband of the width A, e(p) = — (A/2) cos(pa/h), 
we arrive to the expression e(t) — ^'kLo ^ 2k cos(2ku;t), 
where Cik = —AJ2k(eaEo/huj) for k > (J n {x) are 
the Bessel functions). It shows that the electron energy 
within the miniband varies with frequencies which are 
some even harmonics of to: uj^ ven = su> (s = 2, 4, 6 . . .). 
If a bias Ed C is also included to the pump field, e(t) os- 
cillates with two combinations of frequencies ojb ± u>^ ven 
and lub ± io° dd , where ui° dd = slo with s = 1,3,5.... 
However, in the presence of collisions the oscillations 
with Bloch frequency decay, whereas energy oscillations 
with the frequencies imposed by ac field survive. The ef- 
fective electron mass in the nonparabolic miniband also 
varies periodically with the frequency of energy oscilla- 
tions. Now let us suppose that additionally a weak ac 



field E pr — E\ cos(wii + </>) is also applied. The fre- 
quency of this probe field lj% is fixed by an external cir- 
cuit (resonant cavity). Since electron transport in the 
band depends on an instant value of the effective elec- 
tron mass, one should expect the parametric resonance 

for U>e = 

. .even ~~ , .odd 
iO = 01 CO 



lu>\ (I is an integer and lo^ stands for either 
). The most strong parametric resonance 
occurs when I = 2, that is for Loi^ /2 = u)\ @. As in other 
parametric devices |, the parametric resonance due to 
Bragg reflections can result in a regenerative amplifica- 
tion of the probe field. However, currents at harmonics 
of the pump ac field are generated in SL due to strong 
nonparabolicity of its miniband If the parametric 
amplification arises at the same frequencies as the fre- 
quencies of generated harmonics, the effect of harmonics 
blurs out the weaker (oc E%) effect of small-signal gain. 
This problem is well-known for the parametric ampli- 
fication in Josephson junctions, which also have strong 
nonlinearity 0]. 

We are interested in manifestations of the parametric 
resonance due to Bragg reflections in the presence of col- 
lisions, i.e. in the miniband transport regime [5j. Here 
two main questions arise: Can the parametric resonance 
provide a high-frequency gain in the miniband transport 
regime? Is it possible to avoid space-charge instability? 
Some of these problems have been discussed earlier. In 
1977 Pavlovich first used Boltzmann transport approach 
to calculate the coefficient of intraband absorption of a 
weak probe field (o>i) in SL subjected to a strong ac pump 
of commensurate frequency (w) 0] . He briefly mentioned 
a possibility of negative absorption for some u>i/u>. How- 
ever, neither physical origin of the effect nor its compati- 
bility with conditions of electric stability were addressed 
in this pioneer work. Further, in a recent letter (Io| . we 
presented numerical support for a possibility of paramet- 
ric amplification without formation of electric domains in 
the miniband transport regime. Solving numerically bal- 
ance equations for SL [ll[ we demonstrated a feasibility 
of gain at even harmonics. In this situation, we observed 
that gain can exist in the absence of NDC. It guarantees 
electric stability for moderate concentrations of electrons 

m. 
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FIG. 1: (color online) Two schemes of the parametric am- 
plification in superlattice without corruption from generated 
harmonics. In the presence of ac pump (red online) of the 
frequency u), parametric gain for a weak signal (green online) 
of the frequency uj\ arises either at uj\ =: w/2, 3w/2, . . . in 
biased SL (Edc 7^ 0) or at uii =: 2lu,Alo, ... in unbiased SL 
{E dc = 0). 

In this paper, we analytically calculate gain of a weak 
high-frequency (THz) probe field in SL miniband under 
the conditions of parametric resonance, Ue /2 = u>i, 
caused by the action of a strong ac pump field. The 
physical origin of the parametric resonance is a periodic 
variation of effective electron masses in miniband and, at 
high THz frequencies, also a variation of specific quan- 
tum inductance. We prove that for a proper choice of 
relative phase <p a power is always transferred from the 
pump to the probe field. Furthermore, we show that the 
same pump field also modifies free carrier absorption in 
SL. We find that the gain caused by the parametric reso- 
nance can sufficiently overcome the modified free carrier 
absorption and simultaneously remain unaffected by the 
generated harmonics of the pump only in two distinct 
cases: For amplification at half-harmonics in biased SL 
and for amplification at even harmonics in unbiased SL 
(Fig. [I} . In both these cases we predict a significant am- 
plification at room temperature in the absence of NDC. 

Within the semiclassical approach |5| we first solved 
Boltzmann transport equation for a single miniband and 
bichromatic field E p (t) + E pr (t) with commensurate fre- 
quencies. Then we calculated the phase-dependent ab- 
sorption of the probe field, which is defined as 

A=(V(t)cos(u 1 t + 4>)) t , (1) 

where V(t) — V(t)/V p is the electron velocity de(p)/dp 
averaged over a distribution function satisfying the Boltz- 
mann equation and (. . .)t means averaging over a time pe- 
riod which is common for both pump (ui) and probe (u>i) 
fields. Gain corresponds to A < 0. Note that through 
the paper the averaged velocity V, averaged energy W — 
e/\e eq \ and field strengths £0,1 and Ed c are scaled to the 
Esaki-Tsu peak velocity V p — (Aa/iH)fi(T), the equilib- 
rium energy in absence of fields e ea = — (A/2)/x(T) [5| 
and the critical field E c = H/ear [2|, respectively. The 
temperature factor is fi(T) = h(A/2k B T) / I (Aj2k B T) 
(here Io,i(x) are the modified Bessel functions) [l3j]. Ab- 



sorption of a weak (E\ — > 0) probe field in SL (Eq. [T]) is 
linear in E\. It can be naturally represented as the sum 
of phase-dependent coherent and phase-independent in- 
coherent components A = A co h + A inc . 

Parametric effects in the absorption are described by 
its coherent component. It has the form 

A coh = -{(3 1 /A)Bcos[2(<t>~4> op t)l (2) 

where the amplitude of coherent absorption B > and 
0i = Ei/{oj\t). The coherent component always pro- 
vides gain if \<f> — 4> op t \ < 7r/4. Gain has maximum at an 
optimal phase 4> op t . Under the action of pump field, such 
energy storage parameters of SL as the energy of elec- 
trons in miniband W and mesoscopic electric reactance, 
which is described by the reactive current I sln oc V sln , 
are simultaneously harmonically modulated. The vari- 
ables B and (j) op t can be represented in terms of the spe- 
cific harmonics of W(t) and out-of-phase component of 
electron velocity V sln (t) as 

B = [Bf f + Bl f ] 1/2 , tan(20 opt ) = -B hf /B lf , (3) 
Bi f = 2Wr(u B ), 

B hf = W s co >s + wi) - 2VF s co >s) + W s co >b - ui) 

- [vr (wb+wi)- vr {u B - wo] , (4) 

where the index s is the same as involved in the condition 
of parametric resonance and the Fourier components of 
the quantum reactive parameters are given by 

w r = -J2 J iWi J i-^) + J i+k(P)}K(oj B +iu>), 
1 

V k n = Y, J '(P)[Ji+kW)-Ji-k(f3)]K(u; B + lu J ), (5) 
1 

w t = - Y. j W)I j i-^P)- j i+M vET ^b + iu). 

In Eqs. [U P = E /(wt), the Esaki-Tsu drift veloc- 
ity V et {lu b ) = T+j&E @' ^ and Esaki " Tsu enef gy 
K(u>b) — i + (J" BT )ia [a 111! determine the dependence of 
W£ os (u B ), wfi(uj B ) and V£ in (w B ) on the dc bias E dc . 
It worth to notice that instead of harmonics of energy we 
alternatively can consider harmonics of effective electron 
mass because m _1 (e) oc W. 

In the low frequency range ujt, wit < 1, we found that 
Bhf — > and therefore B = Bif, while for THz fre- 
quencies (ujt > 1) both terms Bif and Bhf contribute to 
B. The behavior of the absorbtion amplitude B at THz 
frequencies has two peculiarities. First, influence of the 
out-of-phase component of electron velocity at the pump 
frequency and its harmonics also becomes important. As 
follows from Eq. ([5|) , it describes inductive response of in- 
ertial miniband electrons to ac field in the limit lot ^> 1: 
= E q /utL, L- 1 = 2J Q (l3)J 1 (l3)p- 1 K{u B ) [3. 



Second, interaction of miniband electrons with THz fields 
has quantum nature [H}. Therefore, even a very weak 
probe field produces a back action on the SL reactive 
parameters. This is indicated by an appearance of vir- 
tual processes of absorption and emission of one quan- 
tum of the probe field (±fiu)i) in the expression for Bhf 
(Eq. Q). In particular, Bhf is determined by the dif- 
ference between changes in electron energy at absorption 
W(ujb+ui) — W(lob) and emission W(lob) — W(ujb— Wi). 
The asymmetry in the elementary acts of emission and 
absorption is caused by scattering. It resembles corre- 
sponding asymmetry revealed in the quantum description 
of THz Bloch gain in dc biased SLs [15J] . 

We turn now to the analysis of the incoherent compo- 
nent of absorption Aj nc , which is independent on both 
the ratio uj\/uj and phase difference <f>. It can be repre- 
sented as 

A vnc = y [Vdc(UB + Wl) - Vdc(uB - Wi)] , (6) 

where Vdc = (V)t is the drift velocity induced in SL by 
the pum p fi eld alone. It is determined by the well-known 
formula [16[ 



jfiP^'iLUB+lio). 



(7) 



A.inc describes the free carrier absorption modified by the 
pump. Naturally, A;„ c becomes the usual free carrier ab- 
sorption Ai nc ex (I+o-^t 2 ) -1 in the absence of pump field 
(i?o = Edc = 0). Remarkably, as follows from Eq. ([6]), 
the pump could suppress the free carrier absorption (if 
Vdd^B + w i) ~ Vdd^B — w i)) or even make its value 
negative (if V dc {w B ~ u>i) > V dc (u B + wi)). 

On the other hand, it is easy to see that in the qua- 
sistatic limit, lj%t <C 1, the finite difference in Eq. (J6]) 
goes to the derivative dV dc /dE dc , which determines the 
slope of dependence of V dc on dc bias at the working 
point Edc- The sign of this derivative controls electric 
stab ility against spatial perturbations of charge density 
[l2l 17 1 : For negative slope dV dc /dE dc < destruc- 
tive space-charge instability arises inside SL. In contrast, 
dVdc/dEdc > is the necessary condition for absence of 
the electric domains in moderately doped SLs [l2| ■ 

For general case loit > 1 our numerical analysis 
showed that the sign of finite difference (JHJ) is almost al- 
ways same as the sign of the derivative dV dc /dEd c , if SL 
is unbiased (Edc = 0) or only weakly biased. Therefore 
Ainc > guarantees electric stability. The total absorp- 
tion A = A co h + Ai n c still can be negative in conditions 
of electric stability if \<j> — 4> op t \ < 7r/4 and | A co h |> A inc - 
In the case of unbiased SL (Fig. [1]), such situation is il- 
lustrated in Figs 12131 Fig. [2] shows the regions of negative 
absorption (A < 0) at even harmonics together with the 
regions of NDC (dVd c /dEd c < 0) in uE plane. Here 
the phase is chosen to be optimal (Eq. [3]). The values 
of Eq and u> resulting in electric instability (red areas in 
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FIG. 2: (color online) Amplification at even harmonics in 
unbiased superlattice for cj> = <j>o P t- Regions above the marked 
curves correspond to gain at ui =: 2ui, \lo, 6uj. Dark (red 
online) areas correspond to electric instability. 
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FIG. 3: (color online) Magnitude of negative absorption at 
even harmonics (marked curves) as a function of the pump 
frequency ui for the fixed pump amplitude Eq — 5.1 and (j> = 
(f> op t- Dark (red online) segments indicate intervals of NDC. 



Fig. [21) are close to the lines of Bessel roots Jq(/3) = 0. It 
can be explained noticing that for Edc - > transition to 
NDC is accompanied by absolute negative conductivity 
(ANC) [l2j. However, as can be derived from Eq. [71 in 
the limit E dc -> 0, ANC arises only for J (j3) ~ HU. 
Importantly, the regions of gain and areas of instability 
overlap only in limited ranges of the pump amplitudes 
and frequencies. Moreover, the magnitude of domainless 
gain is significant even at room temperature (Fig. [3]) . To 
estimate gain a in units cm -1 [l5l ] we used the formula 
a = a x (A/ Ex) with a = 8neNVp/ (E c n r c) and the fol- 
lowing typical semiconductor SL parameters: a = 6 nm, 
A = 60 meV, electron density N = 10 16 cm" 3 , r = 200 
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waveguide arrays [l8| and microcavity SLs [19] for light, 
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FIG. 4: (color online) Amplification with a low threshold at 
uji = lo/2 in biased superlattice for Edc = 1 and (f> — (f> op t- 
Marked region between curves corresponds to gain, while dark 
(red online) area corresponds to electric instability. Inset: 
Magnitude of negative absorption as function of the pump 
amplitude Eq for lot = 0.25. 



phononic microcavity arrays 
in perpendicular electric field 
lattices for ultracold atoms [22 1 



20], carbon nanotube SLs 
2l| , and dissipative optical 
These SLs were specially 



suggested and designed to manifest effects of Bloch os- 



cillations hj, Ml 21 



|22j or ac field [H LIU in a sin S le 
band and therefore potentially can be used to observe the 
parametric amplification. 

In summary, we described physical mechanisms for the 
parametric resonance and resulting high-frequency am- 
plification in an energy band. The parametric amplifica- 
tion of a weak signal is possible without negative differ- 
ential conductance. Parametric effects due to Bragg re- 
flections in ac-driven lattices are no less important than 
manifestations of Bloch oscillations in the case of a pure 
dc bias. 

This work was supported by Academy of Finland, Emil 
Aaltonen Foundation, and AQDJJ Programme of ESF. 



fs, refractive index n r = V13 (GaAs) and T — 300 K. 

Even harmonics of the pump satisfy the parametric res- 
onance condition Luf en /2 = u>i. For unbiased case only 
this scheme provides amplification which is unaffected 
by generated harmonics. On the other hand, for E c i c =/= 
subharmonics of the pump (lo\ =: w/2, 3w/2, . . .) sat- 
isfy another parametric resonance condition u>° dd /2 = u>i. 
We found that regions of gain at different half-harmonics 
and areas of electric instability (NDC) have no overlap- 
ping for many values of Eq and lot. Fig. 0] illustrates 
this for amplification at lo\ — lo/2 and E^ c — 1. Here 
threshold is very low while gain is still significant even 
at Eq < 0.5 (Fig. [4l inset). We explain it analyzing the 
behavior of both A co h and A{ nc for small Eq. First, for 
lui/lo = 1/2 relatively large first harmonics (s = 1) of 
the quantum reactive parameters contribute to A co h < 
(Eq. 2]). Second, the tangent to the curve describing a 
dependence of Vdc on Ed c (Eq. [7]) has a small positive 
slope at the working point Ed c — 1- Following Eq. [6] it 
results in a rather small Ai nc > 0. Therefore the total 
gain A < is not small. 

In this Letter, we focused on the phase-sensitive degen- 
erate parametric amplification of THz fields in SLs. Our 
theory can be directly extended to describe nondegener- 
ate phase-insensitive amplification. Here, at least for the 
case of unbiased SL, regions of NDC in loEq plane are still 
located only near Bessel roots lines (c/. Fig. [5J). There- 
fore, by a proper choice of amplitude and frequency of the 
pump it is also possible to reach electrically stable ampli- 
fication of weak signal (loi) and idler (L02) fields satisfying 
the parametric resonance condition uj^ ven — lo\ + lo2- 

The parametric effects in a nonparabolic energy band 
should exist not only in semiconductor SLs but also 
in other artificial periodic structures, including periodic 
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